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Abstract 

We study a stochastically-driven standard map. The addition of a noise term destroys the 
invariant manifolds that organize the phase space which allows for more widespread transport 
than in the noiseless case. Using appropriately defined hitting times to quantify the dynamics, 
we identify two qualitatively different classes of transport: linear and nonlinear. Linear transport 
is primarily driven by the stochasticity in the system, while nonlinear transport results from a 
combination of the nonlinear dynamics and stochasticity and provides a significant speed-up in 
transport. 


1 Introduction 


Symplectic maps represent the discrete-time analogues of Hamiltonian dynamical systems [Mei92] 
which have a wide range of applications, for instance in the study of magnetic fields |Mor00] and 
fluid dynamics |Are84] . Such systems often display rich dynamical phenomena including chaotic tra¬ 
jectories, periodic orbits, and invariant circles, making them a popular subject of interdisciplinary 
research for mathematicians, physicists, and other scientists [AA90] . A particularly well-studied 
example is the two-dimensional area-preserving map on the cylinder T x M given by 


/o : 


x' = X + y' mod 1, 

y' = y + (^9{x). 


( 1 ) 


Here x and y represent angle and action variables, respectively, x' and y' denote their values at 
the next time iterate, e is a nonlinearity parameter, and the force g is assumed to be smooth and 
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periodic (mod 1). In particular, the choice 


9{x) 


— sin(27rx) 

ZTT 


(2) 


yields Chirikov’s standard map |Chir79] . which exhibits a rich ensemble of dynamical phenomena 
despite its simple form and has therefore been the subject of significant research. 

When e = 0 the phase space of the standard map, Eqs. m and (121) . is foliated by invariant circles 
with constant action on which the dynamics is a simple rotation with rotation number lo = yo- KAM 
Theory guarantees the persistence of invariant circles in the standard map under small perturbation 
when u! is “sufficiently irrational” [dinu] . For moderate values of e the dynamics of the standard 
map can be characterized as either resonant or non resonant. In the resonant regions, those around 
where p-o; = q, {p, q) G Z^, contractible circles, often referred to as secondary circles or islands, arise 
alongside chaotic orbits |Dua08] . These effects are most pronounced near the low-order resonances, 
or those with small {p, q). The dynamics of the non resonant regions is similar to the dynamics 
when e = 0. Invariant rotational circles, homotopic to the circles in the e = 0 map, permeate this 
space, however thin secondary circles and chaotic orbits do exist between them. 

The invariant circles of Eqs. m-m, as well as many other systems, organize the phase space 
and determine the extent of possible transport [MMP841 |SCL-|-09l . Each invariant circle acts as a 
barrier for motion - insulating dynamics on either side from one another |GMS86j . Thus, in the 
noiseless case transport is limited to the movement of trajectories throughout a single invariant 
set. More widespread transport can only occur once one or more circles acting as barriers are 
destroyed as a result of an increase in the nonlinearity parameter e |BK84| . Once destroyed, a circle 
typically gives rise to a cantorus - a fractional dimensional cantor set - that allows a slow “leaky” 
transport [BM931 rLB86[ IMM921 IPerSO] . The breakup of these invariant circles in the standard map 
and other area- and volume-preserving maps has been an active area of research |Gre79l IKM891 
IMac92[lMP^IFM14] . 


However, all real-world systems - either natural or man-made - display some degree of noise or 
stochasticity, usually in the form of some temporal fluctuations. In many classes of physical systems 
it has been shown that even very small amounts of noise can fundamentally alter the dynamics of 
the svstem [BG961IZK02]. Investigation into the effect of stochasticitv in the context of Eos. © and 
other well-studied area- and volume-preserving maps to date have been limited |Fro75l IKRW82] . In 
this paper we study the effects of an added stochastic term in Chirikov’s standard map. Specifically, 
we assume the sinusoidal form of the force g, and, after adding a stochastic term we obtain the 
stochastic standard map 


fa-- 


x' = X + y' mod 1, 

y' = y + sin(27rx) 
Ztt 


(3) 


Here, is an iid random variable generated at each time iterate. For simplicity, we assume here 
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that is drawn from the zero-mean normal distribntion with variance cr^, a tunable parameter 
that quantifies the noise intensity. We note that the particular shape of the noise distribution 
does not alter the results we present below, provided that the variance of the noise is but for 
simplicity we consider Gaussian white noise, i.e., ~ 

The addition of the stochastic forcing term dramatically alters the dynamics. Most notably, 
the invariant sets that are so pivotal - fixed points, periodic orbits, and circles - are broken. An 
immediate consequence is that the action, i.e., the y variable, of all orbits will be unbounded. It 
is our goal to understand how transport occurs in this new paradigm, and therefore we will study 
the dynamics in the action (or y) direction. 

In this paper we employ appropriately defined hitting times to quantify the behavior of the 
stochastic standard map ([3|) and find that the transport is a novel combination of the linear noise 
and the dynamical nonlinearity. In particular, when the nonlinearity (i.e., the parameter e) is 
small, we show that the transport is dominated by the noise term and can be well-captured by 
simple Brownian motion properties. However, for larger nonlinearity we observe that the noise 
combines with the dynamics to give rise to regions of rapid nonlinear transport. Specifically, there 
is a significant speed-up in transport in the resonant regions of phase space relative to the non 
resonant regions. 

The remainder of this paper is organized as follows. In Sec. [2] we present a brief survey that 
demonstrates the effect of temporal noise on the dynamics and define the hitting times we use for 
quantifying the transport dynamics. In Sec. [3] we begin our analysis, first presenting numerical 
results for hitting times, then presenting an approximation using a simple Brownian motion. In 
Sec 0] we illustrate the scaling properties of the transport and highlight the nonlinear transport 
effects in the system. In Sec. [S]we present an additional demonstration of the nonlinear transport 
effects. Finally, in Sec. [6] we conclude with a discussion of our results. 

2 Survey and Definitions 

We begin by demonstrating the effect that added stochasticity has on the dynamics of ([3|). In 
particular, we consider three levels of noise: a = 0 (i.e., no noise), 10“^, and 10“^. Setting the 
nonlinearity parameter to e = 0.5 we simulate ([3|), using several different initial conditions for 
each value of a and plot the results in Fig. [TJ Results for <7 = 0, 10“^, and 10“^ are shown in 
panels (a)-(c), respectively, and different colors indicate different orbits obtained using different 
initial conditions. Each orbit plotted consists of 800 time iterations. The a = 0 case [panel (a)] 
corresponds to the classical standard map, i.e., Eqs. dB)-® where each orbit shown is bounded 
away from every other. When cr is increased to 10“^ [panel (b)] we observe a similar structure 
to the phase space, save for a slight thickening of the orbits. While mixing will eventually occur, 
the short-term dynamics shown in panel (b) resemble to a remarkable degree those in panel (a). 
When a is further increased to 10“^ the thickening of each orbit is made even more pronounced. 
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While the orbits are not as easy to distinguish as in panels (a) or (b), the overall structure of the 
dynamics is maintained to a certain degree. 



Figure 1: Effect of noise. Phase space of the stochastic standard map ([3]) for 
fixed e = 0.5 and noise intensities cr = 0 (a), lO”"* (b), and 10“^ (c). 


This survey clearly illustrates the fundamental change in dynamical behavior that occurs in the 
standard map with the addition of stochasticity. In particular, the invariance that is key to the 
dynamics in the noise-less standard map is destroyed. Most importantly, this allows for transport 
not just within an invariant set in the noise-less case, but across the entire phase space. Evidence for 
this more robust transport can be first observed in panel (b) (cr = 10“^) with the slight thickening 
of the orbits. However, transport is more clearly seen in panel (c) {a = 10“^) as orbits are already 
overlapping after only 800 iterations. 

To investigate the effect of stochasticity on the new dynamics of ([3|) and especially to quantify 
the transport throughout the phase space, we will use an appropriately defined hitting time. In 
particular, we examine how long it takes an orbit of the stochastic standard map to reach a distance 
a from an invariant set of the noiseless standard map. Our exploration of these dynamics will be 
split into two parts. In Secs [3] and [3] we examine the behavior near a rotational circle of the noiseless 
system and in Sec [5] we study transport near a periodic orbit of the noiseless system. 

More precisely, let 


^o{xo,yo) = {fo{xo,yo)\t €N}. (4) 

denote an orbit of the noise-less map /o [see Eqs. ©]• We examine the distance between <l>o(a;o,yo) 
and orbits of the stochastically forced standard map ([3]) with nose level a beginning at the same 
initial condition (xo,yo)- Specifically, we calculate for a given distance a the hitting time r^, which 
is defined as the first time t that the orbit of fa [see Eqs. ([3])] beginning at (xq, yo) equals or exceeds 
a distance a from the set 4>o(xo,yo): 

Ta{xo, yo] cr) = inf{t G N | d[/* (xq, yo), ^o{xo,yo)] > a}- (5) 
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Thus, given an appropriate distance metric d{-, ■) on the space T x M, the hitting time represents 
the time it takes for a noisy trajectory to diverge from the noiseless trajectory by a distance a, and 
here we are interested in the expected hitting time. 

As an illustrative example to begin with, in Sec. [3] we draw initial conditions from the golden- 
mean circle, the rotational circle with ui = , which we illustrate in Figs. [2]^ a), (b), and (c) for 

nonlinearity parameters e = 0.01, 0.2, and 0.5, respectively. The golden-mean circle if denoted as s 
thick blue curve, and other orbits are denoted in red. We choose the golden-mean circle in particular 
because of its robustness, i.e., the golden mean circle is believed to be the last circle to survive and is 
only destroyed when the nonlinearity parameter is increased to e = Cc ~ 0.971635 [Mac93[ IOP08] . 
Points on this circle are computed using the quasi-Newton method developed by de la Llave et 
al. [HdlLSMllFMTi] . 



Figure 2: Golden-mean circle. For the noise-less standard map, the golden- 
mean invariant torus for e = 0.01 (a), 0.2 (b), and 0.5 (c) as the thick blue 
curve, compared to other trajectories in red. 


Given the topology of the cylinder and the fact that we are most interested in exploring the 
noise-induced transport in the action variable, we consider the simple distance metric defined solely 
by the displacement in the y-direction. In practice we begin with an initial condition on the circle 
and iterate /o 1000 times to generate a sample of points on d>o. To measure the distance from a 
point (x, y) to <ho we use a linear interpolation to approximate the point on the circle {xc, Vc) with 
the same x coordinate, Xc = x. The distance is then given by d = \yc — y\- This is well-defined 
since every rotational invariant circle of the standard map is a graph |Mei92| . This choice for d is 
further motivated by a simplification that it will allow in the analysis we present below. 


3 Hitting Times and Brownian Motion 

Using the hitting times defined above we now study the effect of different noise levels on the 
dynamics of ([3]). Starting on the golden mean circle described above, we calculate the hitting 


5 







































times Ta for various distances a, nonlinearity parameters e, and noise levels a. Since the process 
is stochastic, we will be interested in the expected hitting times, and therefore we will calculate 
for each set of parameters the mean Ta from 10^ realizations. In Fig. [3] we plot the hitting times 
Ta VS the distance a for e = 0.01, 0.2, and 0.5 in panels (a)-(c), respectively, and for noise levels 
a = 10“^, 10“^'^, and 10“^, plotted in blue, red, and green circles, respectively. For each value 
of e, as expected, the expected hitting time Ta increases both as the distance o increases and as 
the noise level decreases. Note, however, that for smaller e, e.g., panel (a), the rate at which Tq 
increases with a is quite regular, increasing very much like a power-law, while for larger e, e.g., 
panel (c), the rate of increase is much less regular, especially at larger distances a. 
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Figure 3: Hitting times. For nonlinearity parameters e = 0.01 (a), 0.2 (b), 
and 0.5 (c), we plot the mean hitting time vs distance a for noise intensities 
a = 10“^, 10“^-®, and 10“^, which are plotting in blue, red, and green cir¬ 
cles, respectively. Each data point represents the mean over 10^ realizations. 
The Brownian motion approximations Ta = c? jcP' are plotted as appropriately 
colored dashed curves. 


To gain a better understanding of the behavior of we consider the behavior of the process 
in the limit of small nonlinearity, i.e., e —>■ 0. In particular, we note that in this limit the vertical 
motion depends only on the noise terms, i.e., y' — y = Thus, given our notion of distance defined 
by displacement in the action variable, we ignore motion in the angle variable, obtaining effectively 
a one-dimensional system. Next, since the noise is assumed to have variance a simple 
application of the Central Limit Theorem implies that, at a large enough time t, the distribution 
for the displacement yt — yo is Gaussian with variance cr^f. It follows that at large enough times 
the discrete process can thus be approximated as finite-time slices of a one-dimensional Brownian 
motion with variance for which the hitting time is well known to be given by 


Ta = 



( 6 ) 


Eq. Q thus provides an approximation for the hitting times for small e, assuming that is not 
too small (i.e., a and a are not simultaneously too big and small, respectively). 
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To compare this approximation with our results, we plot in Figs. [3)^a)-(c) the Brownian motion 
approximation Ta = as dashed curves. For the case of e = 0.01 [panel (a)], the approximation 

captures the observed behavior extremely well as the scaling Ta oc is almost exactly obtained. 
For e = 0.2 the approximation remains remarkably accurate for a not too large, until it begins to 
break down near a ~ 0.4. Finally, for e = 0.5 the approximation breaks down earlier still, near 
a ~ 0.05, but provides a good benchmark for smaller a. We note that in our simulations there is a 
small discrepancy for a = 10^ and small a, as the green circles towards the beginning of each plot 
slightly overshoot the approximation. We find that this is a result of modeling a discrete process 
with a continuous one - an effect that arises when the expected hitting time is not large enough. 
These results beg the question of why the approximation fails at certain distances a for larger 
nonlinearity parameters e - a point we address next. 

4 Rescaling and Nonlinear Transport 

While the approximation via Brownian motion presented above represents a useful benchmark 
for understanding the dynamics of the noisy system, two interesting points remain. First, as 
the Brownian motion approximation fails, we observe a speed-up in the transport of the system. 
In other words, when the approximation Tq ~ loses accuracy, the observed mean hitting 

time is always smaller, indicating that transport always occurs quicker - and never slower - than 
predicted by the Brownian motion approximation. Second, the speed-up in transport of the system 
seems to scale with the noise level. In particular, for a given value of e, the deviation from the 
approximation appears the same up to a rescaling for different values of a. [See in particular the 
right-hand-side of Fig. Ha), where the blue, red, and green circles appear to undershoot their 
respective approximations by the same amount in each case.] 

We begin by noting that the dependence of the hitting time on noise level can be scaled out of 
the Brownian motion approximation in Eq. (I6|) by considering the rescaling Ta •—> cP'Ta- Therefore, 
we plot in Fig. 0] the scaled quantities cr^Ta vs the distance a using the same results presented 
in Fig. [3l Results using e = 0.01, 0.2, and 0.5 are plotted in panels (a)-(c), respectively, and 
results using o = 10“^, 10“^'®, and 10“^ are plotted in blue, red, and green circles, respectively. 
As expected, for small e, e.g., panel (a), the results collapse onto the curve o^Ta = c?■ More 
surprisingly, however, we observe that the results for larger e, i.e., panels (b) and (c), also collapse 
on one another, even in the regime where the Brownian motion approximation fails. 

The fact that even the deviations from the Brownian motion approximation scale is remarkable 
and suggests two important points. First, the observed speed-up in transport is not solely a result of 
the added stochasticity, but also from the underlying nonlinear dynamics which are deterministic. 
Second, the transport we observe falls into two regimes. We refer to these two regimes as a linear 
transport regime and a nonlinear transport regime. The linear transport regime corresponds to 
dynamics that adhere well to the Brownian motion approximation, indicating that the dynamics 
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Figure 4: Scaled hitting times. For nonlinearity parameters e = 0.01 (a), 0.2 
(b), and 0.5 (c), we plot the scaled mean hitting time a^Ta vs distance a for 
noise intensities a = 10“"^, 10“^ ®, and 10“^, which are plotting in blue, red, 
and green circles, respectively. Each data point represents the mean over 10^ 
realizations. The Brownian motion approximation cr^Ta = is plotted as a 
black dashed curve. 


are driven primarily by the linear stochastic term and subsequently the expected hitting time is 
well-described by the power-law Tq oc a^. The nonlinear transport regime corresponds to dynamics 
that fail to be described well by the Brownian motion approximation and where the underlying 
nonlinear dynamics of the standard map contribute to the significant speed-up in transport. 

Finally, these results beg for the answer to the question of what the underlying cause for the 
transition from linear to nonlinear transport. While we find that this effect depends significantly 
on the structure of the dynamics, and thus on both the nonlinearity parameter and the initial 
conditions, we find that transport is linear in regions of phase space far from low-order resonance 
which are dominated by rotational circles. Conversely, transport is nonlinear in resonant regions of 
the phase space which are dominated by secondary circles. This can be understood as follows. In 
regions of phase space primarily populated by rotational circles the dynamics push the trajectories 
primarily in the angle direction, so that the majority of the motion in the action direction is driven 
by the stochasticity and is captured well by the Brownian motion approximation. On the other 
hand, in regions of phase space primarily populated by secondary circles, the underlying dynamics 
can induce a significant displacement in the action variable in just a few iterations. Thus, in these 
regions transport is facilitated by the nonlinear dynamics which can transport much quicker than 
the stochasticity can. 

As an example, consider the case of e = 0.5, the hitting times and scaled hitting times for which 
are presented in Figs. He) and He) • We observe that the deviation from the Brownian motion 
approximation, i.e., the kink in the results, occurs roughly at o ~ 0.05. In Fig. [2] we can see that 
this is approximately the distance from the golden mean circle to the resonance centered around 
y = 0.5. Another sizable kink occurs at a ~ 0.25, which corresponds to the distance from the 
golden mean circle to the resonance centered around y = 1. 
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5 Nonlinear Transport: Another Example 


To further demonstrate the nonlinear transport effects in the dynamics of Eqs. ([3|) we consider 
a modified scenario where, rather than starting at the golden mean circle or another rotational 
circle, we begin at the period-two orbit ^>0 = {(0, ^), (^, ^)}. This orbit is plotted as blue dots in 
Figure [5l a), with other trajectories plotted in red, for e = 0.2. Conveniently, this orbit is fixed for 
every value of e, making it an illustrative choice of initial conditions for the following experiment. 



e = 0.2 e = 0.2 



Figure 5: Example two. (a) The orbit $o = {(0, 5)1 ^)} plotted in blue dots, 

compared to other trajectories in red for e = 0.2. (b) Scaled hitting times a^Ta 
vs distance a for noise intensities cr = 10 “^, 10 “^’®, and 10 “^, plotted in blue, 
red, and green circles, respectively. Each data point represents the mean over 
10^ realizations. The Brownian motion approximation ct^Tq = is plotted as 
a dashed black curve, (c) A zoomed-in view on the flat region in panel (b) with 
finer resolution. 


We proceed by simulating the dynamics of Equations ([3]), again calculating the mean hitting 
times Ta as a function of the distance a. We note that, since the orbit is not a graph over the angle 
variable, we calculate distance simply as the y-displacement from the value yo = 0.5. In Figure[5jb) 
we plot the mean hitting times scaled by the noise squared, a'^Ta vs the distance a for cr = 10“^ 
(blue circles), 10“^'^ (red circles), and 10“^ (green circles), again averaged over 10^ realizations. 
The Brownian motion approximation a'^Ta = is denoted by the dashed black curve. Note that 
the results for different values of a collapse nicely when scaled appropriately. 

In Figure EKb) we observe two nonlinear transport effects, manifesting as significant deviations 
from the linear transport along the Brownian motion approximation. First, and most subtle, we 
observe that for small distances (i.e., a < 10“^) there is a speed-up in the hitting times. Upon 
further investigation, we find that this turns out to be due to the curvature of the trajectories 
around the period-two orbit depicted in Figure [5Ka). In particular, the mildly sinuous shape of the 
orbits - an effect of the nonlinear dynamics of the system - causes a small speed-up in transport 
with respect to the action variable. Second, we see a more pronounced nonlinear transport effect 
at larger distances (i.e., the second- and third-to-last data points) where the hitting times flatten 
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out. Since the logarithmic scaling of the plot may diminish the apparent size of this region, we plot 
in Figure El^c) a zoomed-in view with a ten-fold finer discretization. Here we see more clearly the 
dramatic, flat section of the plot in the range 0.4 < a < 0.6 where transport occurs very quickly. 
Looking back to the foliation of the phase-space in Figure E^a), we note that this flat region 
corresponds almost precisely with the distances from the initial value yo = 5 to and extending 
through the large resonances centered at (^,0) and (^,1). Once past these large resonant regions 
the transport returns to the linear regime, as shown by the final point in Figure [5)^b). 

6 Discussion 

In this paper we have studied the dynamics of Chirikov’s standard map with an added stochasticity 
term [see Eqs. ([3])]. The added noise term in the stochastic standard map destroys the invariant 
manifolds present in the noiseless case that organize the phase space and bound transport. The 
destruction of these invariant objects facilitates widespread, unbounded transport in the action 
direction not present in the noiseless case. Using appropriately defined hitting times, we have 
quantified the transport that occurs in the stochastic standard map and found that transport falls 
into two broad categories: linear transport and nonlinear transport. In the case of linear transport, 
movement in the action direction is dominated by the stochastic term and is well-described by a 
simple Brownian motion such that hitting times scale with the square of the distance. In the case 
of nonlinear transport, the stochasticity combines with the underlying nonlinear dynamics of the 
map to facilitate a significant speed-up in the mean hitting times. Importantly, we find that linear 
transport prevails in regions of phase-space dominated by non-resonant dynamics, while nonlinear 
transport prevails in regions of phase-space dominated by resonant dynamics - which become more 
pronounced as the nonlinearity parameter of the dynamics is increased. 

The effect of added stochasticity to the standard map or other conservative systems has been 
studied in a handful of other works, e.g., see Refs. [Fro75[ lKRW82| . however very few recent results 
exist. To our knowledge, this is the first study concerned with the effect that added stochasticity 
has on the the transport that occurs as a result of breaking the invariant objects in phase space. 
Subsequently, our results open new questions for further investigation. First, the examples pre¬ 
sented in this paper were chosen using parameter values such that the noise-less phase space was 
primarily foliated by circles and (periodic) fixed points. One interesting avenue for investigation 
will be to study how transport occurs as larger chaotic seas emerge (i.e., at larger values of the 
nonlinearity parameter). We hypothesize that the emergence of significant chaotic seas will cause 
even more significant speed-up in transport than do resonant circles. Second, we have used as a 
primary example in this paper Chirikov’s standard map due to its simplicity and widespread popu¬ 
larity. However, the effect that noise has on other, possibly more complicated conservative systems 
- discrete or continuous - will be an interesting question to ask. Finally, we find it remarkable that, 
although the invariant manifolds organizing the phase space in the noiseless case are broken with 
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the addition of noise, the dynamics remain remarkably robust. An investigation into the limits of 
this robustness, e.g., in terms of the noise intensity and/or the nonlinearity parameter, would be 
worthwhile. 
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